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ABSTRACT 

ngb-resohitioa  muMichaniicI  seismic  data  collected  from  the  Blake  Rid^  in  the  western 
North  Atlantic  by  the  Naval  Research  Laboratory’s  Deep  Towed  Acoustics/Geophysics  System 
(DTAGS)  show  that  the  bottom  simulating  reflector  (BSR)  in  this  area  is  the  reflection  frtmi  the 
interlhce  between  an  ~440-m-thick  section  of  hydrate-bearing  sediment  overlying  an  — 5-m- 
thick  layer  of  methane  gas-rich  sediment.  The  high  resolution  attainable  by  the  deep-tow 
seismic  system  reveals  normal-fiiuH  oLfrets  of  -20  m  in  the  BSR.  These  growth  faults  may 
provide  a  path  for  vertical  migration  of  methane  initially  concentrated  beneath  the  hydrate- 
bearing  sediment,  enabling  hydrate  to  form  throughout  sediment  above  the  BSR.  Because  the 
BSR  represents  a  methane  gas-methane  hydrate  phase  boundary  rather  than  a  lithologic  or 
diagenetk  horiion,  the  observed  oftet  of  the  BSR  itself  reflects  discontinuities  in  the  pressure- 
temperature  Held  across  the  fault  zones  where  they  intersect  the  BSR. 


INTRODUCTION 

Methane  hydrate,  an  icelike  solid  that 
contains  high  concentrations  of  methane,  is 
found  within  sediments  where  low  temper¬ 
atures  and  high  pressures  are  such  that 
methane  hydrate  is  the  stable  phase  of  a  wa¬ 
ter-methane  mixture.  In  1989,  the  Naval 
Oceanographic  and  Atmospheric  Research 
Laboratory  (now  part  of  the  Naval  Research 
Laboratory  [NRL))  recorded  high-resolu¬ 
tion,  deep-tow,  multichannel  seismic  data 
over  an  area  on  the  Blake  Ridge,  in  the  west¬ 
ern  North  Atlantic,  where  the  presence  of 
methane  hydrate  had  been  inferred  from 
conventional,  surface-tow,  seismic-reflec¬ 
tion  data  and  from  the  drilling  results  from 
Deep  Sea  Drilling  Project  (DSDP)  Sites  102, 
103,  104,  and  533  (Hollister  et  al.,  1972; 
Sheridan  et  al.,  1983).  Evidence  for  the  pres¬ 
ence  of  significant  amounts  of  methane  hy¬ 
drate  from  surface-tow  seismic  data  in¬ 
clude  a  bottom  simulating  reflector  (BSR), 
an  unusually  high  amplitude  reflection  ho¬ 
rizon  that  follows  sea-floor  topography 
(Shipley  et  al.,  1979),  unusually  high  com- 
pressional  velocities  within  the  sediment 
(Dillon  and  Pauli,  1983).  and  anomalous 
Inw-reflectivity  zones  (Dillon  et  al.,  1991). 
These  characteristics  and  evidence  of  nor¬ 
mal  faults  penetrating  the  BSR  arc  shown 
by  NRL’s  high-resolution,  multichannel 
seismic  data  collected  on  the  northwestern 
slope  of  the  Blake  Ridge  (Fig.  1).  The 
Blake  Ridge  is  a  thick  accumulation  of  or¬ 
ganic-rich  Miocene  to  Holocene  terrige¬ 
nous  sediments,  deposited  by  contour  cur¬ 
rents  (MarkI  and  Bryan,  1983),  in  which 
the  upper  -660  m  of  sediment  consists  of 
uniform  hemipelagic  silty  clay  inter¬ 
spersed  with  thin  ( <  I  m)  beds  and  lenses 
of  silty  to  sandy  clay  and  calcareous  mud 
(Hollister  et  al..  1972;  Sheridan  et  al., 
1983). 


HIGH-RESOLUnON  MULTICHANNEL 
SEISMIC  MEASUREMENTS 
The  Deep  Towed  Acoustics/Geophysics 
System  (DTAGS)  is  a  unique  multichannel 


seismic  system  in  which  the  source,  a  Helm¬ 
holtz  transducer  that  generates  a  chirp  from 
250  to  650  Hz,  and  linear  receiving  array  are 
both  towed  300  to  500  m  above  the  sea  floor 
at  full  ocean  depths  (Gettrust  and  Ross, 
1990).  The  array  contains  24  elements  21.0 
m  apart  offset  137-620  m  behind  the  source. 
Shots  are  fired  at  —21  m  (30  s)  intervals.  The 
high-frequency  source  and  close  proximity 
of  the  system  to  the  sea  floor  enable  the  sys¬ 
tem  to  resolve  subbottom  structural  features 
on  a  scale  of  —5  m  vertically  and  —21  m 
horizontally  (Rowe  and  Gettrust,  1993), 
greater  than  the  resolution  possible  from 
surface-tow  multichannel  seismic  systems. 
Normal-moveout-based  semblance  analy¬ 
sis  techniques  were  applied  to  the  common 
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Figure  i.  Map  ahowlng  leeation  of  2.8  km  DTAOS  imiWebannal  aaiamle  Hna  from  Maha  RMga. 
CoMoura  iiNilcata  dapfha  In  malraa;  elrelaa  Indicala  DSDP  sHaa. 
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Figure  2.  Selemle  Mctlon  recorded  on  Bloke  RMge  compiled  from  137  m  offeei  trace  from  each  ahol 
gather.  Bottom  elmulatlng  reflector  (BSR),  -0.63  a  below  eea  floor,  marlia  boae  of  mothane^ydrate 
•tabiiliy  zone.  Two-wey  treveltlme  le  referenced  to  DTAQ8  tow  depth,  -SOO  m  above  eea  floor. 
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shot  gather  data  to  obtain  high-resolution 
(—25  m  vertically,  -  KM)  m  horizontally)  es¬ 
timates  of  the  compressional  vek^ity  within 
the  sediment  as  a  function  of  range  along  the 
ship  track  and  depth  below  the  sea  flixir 
(Rowe  and  Gettrust,  KWS). 

The  seismic  data  discussed  in  this  paper 
were  collected  along  a  2.5-km-long  line  per¬ 
pendicular  to  the  ridge  slope  in  water  depths 
between  3700  and  .3800  m  (Fig.  I ).  The  com¬ 
mon  offset  (137  m)  seismic  section  (Fig.  2) 
shows  numerous  closely  spaced,  flat-lying 
reflection  horizons  caused  by  the  thin  (<4 
m)  layers  of  silty  to  sandy  sediment  and  cal¬ 
careous  clay  observed  in  cores  (Hollister  et 
al.,  1972;  Sheridan  et  al.,  1983).  Below  —0.2 
s  (all  times  given  are  two-way  traveltimes) 
below  the  sea  floor  (bsf),  reflection  ampli¬ 
tudes  decrease  relative  to  the  reflection 
amplitudes  from  the  shallower  reflectors, 
because  of  the  decrease  in  impedance  mis¬ 
match  across  sediment  layers  in  regions 
where  hydrate  has  formed  within  the  sedi¬ 
ment  (Dillon  et  al.,  1991 ).  Normal  faults  dip¬ 
ping  downslope  offset  reflection  horizons 
from  the  sea  floor  through  the  BSR  that  is 
—0.64  s  bsf. 

Compressional  velocity  estimates  ob¬ 
tained  ffom  the  deep-tow  multichannel  data 
(Fig.  3)  show  that  there  is  almost  zero  ve¬ 
locity  gradient  between  the  sea  fltxtr  and 
0.17  s  bsf.  typical  of  shallow,  fine-grained 
marine  sediment  (Gettrust  et  al.,  1988).  Be¬ 
tween  -0.17  and  -0.18  s  bsf.  the  compres¬ 
sional  velocity  increases  sharply,  from 
—  1.65  to  >2.(M)  km/s  (Fig.  3).  Below  this 
depth,  velocity  continues  to  increase,  reach¬ 
ing  2.45-2.70  knvs  at  the  BSR.  These  un¬ 
usually  high  velocities  have  been  measured 
in  hne-grained  marine  sediments  buried 
<6(M)  m  deep  in  areas  in  which  the  fractures 
and  pore  space  are  partially  filled  with  meth¬ 
ane  hydrate  (Mathews  and  von  Huene, 
1985).  Although  the  top  of  this  high-velocity 
gradient  zone  is  at  approximately  the  same 
depth  across  the  section,  comparison  of  Fig¬ 
ure  2  and  3  shows  that  it  does  not  coincide 
with  a  specific  reflection  horizon,  suggesting 
that  the  upper  boundary  of  the  hydrate-bear¬ 
ing  sediment  is  gradational,  as  suggested  by 
Tucholkc  et  al.  (1977).  The  increase  in  com¬ 
pressional  veliKity  with  depth  reflects  the  in¬ 
creasing  percentage  of  sediment  pore  space 
that  contains  methane  hydrate.  Lateral  var¬ 
iability  of  the  compressional  vekKity  within 
the  high-vekx;ity  zone  (Fig.  3)  is  consistent 


FIgurt  3.  Two-dltnvntlonal  Inlarval  valocfty 
modal  darivad  from  OTAGS  data.  Black  area  at 
baaa  of  modal  marks  maximum  depth  for  which 
there  are  reflection  horizons  to  datarmlna  va- 
locltias.  Two-way  trsvaltima  Is  refarancad  to 
DTAQS  tow  depth.  BSR  =  bottom  simulating 
reflector. 
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with  lateral  variability  in  the  sediment  hy¬ 
drate  content. 

The  reflection  horizon  that  is  0.64-0.65  s 
bsf  (Fig.  2).  has  the  high  amplitude  and  in¬ 
verted  polarity  characteristic  of  a  methane 
hydrate-related  BSR  (Shipley  et  al.,  1979). 
This  two-way  traveltime  Is  consistent  with 
that  of  the  BSR  observed  from  surface-tow 
seismic  data  acquired  in  this  area  (Tucholke 
et  al..  1977;  Dillon  and  Pauli,  1983).  The 
deep-tow  seismic  data  reveal  that  this  BSR 
is  a  reflecting  horizon  composed  of  reflec¬ 
tions  from  the  top  and  bottom  of  a  layer  hav¬ 
ing  a  thickness  of  -1-2  wavelengths 
I  Fig.  4A:  Rt>we  and  Gettrust,  1993).  Al¬ 


though  this  layer  is  too  thin  for  estimates  of 
the  compressional  velocity  within  this  layer 
to  be  made  from  the  multichannel  data,  the 
distinct  inversion  in  the  phase  of  the  reflec¬ 
tion  observed  at  the  top  of  the  BSR,  at 
—  1 .29  s  (Fig.  4A),  and  lack  of  conesponding 
density  decrease  (Hollister  et  al.,  1972)  in¬ 
dicate  that  the  compressional  velocity  of  the 
material  immediately  beneath  the  top  of  the 
BSR  is  significantly  lower  than  the  compres¬ 
sional  velocity  of  the  overlying  sediment. 
The  reflection  data  at  the  BSR  were  synthe¬ 
sized  (Fig.  4B)  using  the  full-wavefield 
model  SAFARI  (Schmidt,  1984).  The  best 
fit  to  the  shape  and  amplitude  of  the  BSR 


FIgui*  4.  A:  Part  of  aingla 
traca  from  ahel  gattwr  lo- 
catad  1.0  km  along  tac¬ 
tion.  Normal  phaaa  ra- 
flactiona,  1 . 2.  and  4,  iiava 
paakt  to  right.  Rafloc- 
tlont  with  Invartad 
phaaa,  3,  hava  paaha  to 
loft.  AmptMuda  la  normal- 
Iztd  to  paak  traca  ampli- 
tuda.  B:  Synthetic  aala- 
mogram  that  boat  llta 
traca  In  A  gtnaratod  by 
uaing  comproaalonal  ve¬ 
locity  modal  ahown  In  C 
with  Ricker  wavelet 
aourca.  Relative  ampli¬ 
tude  and  phaaa  of  re¬ 
flected  wavaleta  t-4 
compare  wall  with  data. 
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Amplitude  la  normaUzad  to  paak  trace  ampHtuda.  C:  Beat-Wt  compraaatonal  valeclty  modal  uaad  to 
generate  aynthatlc  aaiamogram  m  B.  Two-way  traveWme  la  ralarancad  to  DTAOS  tow  depth. 
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reflection  was  obtained  using  a  model 
(Fig.  4C)  consisting  of  a  4-m-thick,  low-ve¬ 
locity  (1.3  km/s)  layer  beneath  high-velocity 
(2.4  km/s)  material  and  underlain  by  another 
higher  velocity  (2. 1  km/s)  layer.  Sediment 
compressional  velocities  as  low  as  1 .3  km/s 
may  be  found  in  sediment  in  which  pore 
spaces  contain  free  gas  (Miller  et  al.,  1991). 

Conversion  of  the  seismic  data  from  time 
to  depth  (Fig.  S),  computed  by  using  the  in¬ 
terval  compressional  velocities  derived  from 
these  data  (Fig.  3),  shows  that  the  thin  layer 
causing  the  BSR  reflection  lies  640-660  m 
bsf  and  is  3.S-7.0  m  thick.  The  high  com- 
pressional-velocity  gradient  zone  that  coin¬ 
cides  with  the  top  of  the  zone  of  low  reflec¬ 
tivity,  inferred  to  mark  the  top  of  the 
hydrate-bearing  sediment,  is  -210  m  bsf. 
This  hydrate-bearing  interval  (-210 — 650 
m  bsf)  coincides  with  the  depths  at  which 
high  concentrations  of  methane  were  ob¬ 
served  in  core  samples  from  nearby  DSDP 
sites  (Hollister  et  al.,  1972;  Sheridan  et  al., 
1983). 

FAULT  CONTROLS  ON  METHANE 
DISTRIBUTION 

For  methane  hydrate  to  form  within  this 
— 440-m-thick  layer  of  methane-bearing  sed¬ 
iment,  large  quantities  of  methane  must  be 
available  throughout  the  region.  Below  the 
BSR,  the  thick  sediment  section  contains 
large  amounts  of  rapidly  buried  organic  ma¬ 
terial  (Mark!  and  Biyan,  1983),  a  large  po¬ 
tential  methane  source.  Sufficient  quantities 
of  methane  for  hydrate  formation  may  be 
concentrated  at  the  base  of  the  zone  of  hy¬ 
drate  stability,  marked  by  the  BSR,  by  up¬ 
ward  migration  of  free  methane  gas  or  meth¬ 
ane-rich  pore  fluid,  driven  by  high  pore 
pressures  developed  as  methane  is  evolved 
(Hyndman  and  Davis,  1992).  As  methane 
enters  the  hydrate  stability  zone,  it  forms 
hydrate,  filling  sediment  fractures  and  pore 
spaces  at  the  base  of  this  zone.  The  contin¬ 
uous  BSR  and  unusually  high  compressional 
velocities  throughout  the  sediment  at  the 
BSR  are  evidence  that  sufficient  hydrate  has 
formed  within  the  sediment  to  create  an  im¬ 
permeable  layer  that  inhibits  migration  of 
methane  into  the  overlying  sediments.  Meth¬ 
ane  that  has  continued  to  migrate  from 
deeper  sediment  is  trapped  at  the  BSR.  This 
trapped  gas  may  form  a  reservoir  beneath  the 
impeimeaMe  hydrate-bearing  sediment  layer, 
as  suggested  by  the  full-waveform  model. 

The  normal  faults  resolved  by  these  seis¬ 
mic  data  extend  into  this  methane-rich  sed¬ 
iment  layer  (Fig.  5).  Slope  failure  on  the 
Blake  Ridge,  which  may  be  (Kcurring  along 
this  weak,  gas.sy  sediment  layer  (Dilkm  et 
al.,  l983)or  within  the  non- hydrate-bearing 
sediment  bektw  the  hydrate  stability  zone,  is 
accommodated  in  the  overlying,  more  rigid. 
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hydrate-bearing  sediments  along  these  nor¬ 
mal  faults  (Fig.  5).  The  dip  angle  of  these 
faults,  44°  to  60°,  becoming  shallower  at 
depth,  and  the  increasing  ofeet  with  depth 
are  consistent  with  the  morphology  of 
growth  faults.  Periodic  rupture  along  these 
faults  may  provide  avenues  for  migration  of 
large  quantities  of  methane  past  the  imper¬ 
meable  hydrate-bearing  sediment  near  the 
base  of  the  hydrate  stability  zone  (Reed  et 
al.,  1990).  During  the  time  the  faults  are  in¬ 
active,  hydrate  seals  the  faults  at  the  BSR 
depth,  preventing  methane  gas  from  migrat¬ 
ing.  Movement  along  these  faults  causes  the 
hydrate  within  the  fault  zones  to  dissociate, 
due  to  the  pressure  decrease  and  frictional 
heating  caused  by  the  fault  movement.  At 
this  time,  methane  from  the  reservoir  be¬ 
neath  the  hydrate-bearing  sediment  is  able 
to  migrate  upward  along  the  fault  zone  into 
overlying  non- hydrate-bearing  sediment  ly¬ 
ing  within  the  hydrate  stability  zone.  Re¬ 
peated  fault  rupture  enables  hydrate  to  form 
throughout  a  region  extending  several  hun¬ 
dred  metres  above  the  BSR. 

BSR  FAtLT  OFFSETS 

If  the  BSR  represents  a  pha.se  boundary, 
how  are  the  observed  vertical  offsets  (Fig.  5) 
in  the  BSR  itself  maintained?  It  has  been 
suggested  that  the  BSR  might  be  caused  by 
precipitation  of  carbonate  minerals  in  a  thin 
sediment  layer  at  the  base  of  the  zone  of  hy¬ 
drate  stability  (Lancelot  and  Ewing,  1972). 
However,  the  observed  inverted  phase  of 
the  BSR  reflection  requires  a  negative  im¬ 
pedance  contrast  between  the  sediments 
above  the  BSR  and  the  thin  layer  at  the  BSR. 
A  negative  change  in  impedance  requires  a 
decrease  in  compressional  velocity  or  den¬ 
sity  across  the  interface,  which  is  not  true  of 
the  carbonate  layers  found  at  these  depths 
(Lancelot  and  Ewing,  1972). 

The  abrupt  discontinuities  along  the  BSR, 
coincident  with  the  fault  zones,  may  be 
maintained  by  discontinuities  in  the  temper- 
ature-pres.sure  field  which  exist  across  the 
faults.  Stress  is  higher  on  the  downthrown 
side  of  a  normal  fault  than  on  the  upthrown 
side  (Hafner,  I9SI).  Under  the  higher  pres¬ 
sure.  the  region  in  which  hydrate  is  stable 
would  extend  deeper,  into  warmer  sedi¬ 
ments,  causing  a  downward  shift  in  the 
depth  of  the  BSR  across  the  fault.  The  ob¬ 
served  -20  m  offset  is  equivalent  to  a  pres¬ 
sure  difference  of  ~200  kPa  ( -2  bar)  across 
the  faults,  which  enables  hydrate  to  be  sta¬ 
ble  at  temperatures  - 1  °C  higher.  A  discon¬ 
tinuity  in  the  isotherm  across  the  fault 
caused  by  the  downward  shift  of  colder, 
shallower  sediment  on  the  downthrown 
side,  adjacent  to  the  deeper,  warmer  sedi¬ 
ment,  would  further  enhance  the  apparent 
offiiet.  MacLeod  (1982)  estimated  that  sev- 
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eral  thousand  years  are  required  for  this  iso¬ 
therm  to  regain  equilibrium. 

CONCLUSIONS 

Evidence  from  compressional  velocities, 
reflection  amplitude,  and  reflection  charac¬ 
ter  indicate  that  the  BSR  on  the  Blake  Ridge 
is  the  reflection  from  the  interface  between 
an  ~440-m-thick  hydrate-bearing  sediment 
layer,  which  formed  within  the  zone  of  hy¬ 
drate  stability,  and  an  -4-m-thick  layer  of 
sediment  containing  free  methane  gas  imme¬ 
diately  beneath  it.  Although  hydrate  fotming 
in  the  sediment  at  the  base  of  the  hydrate 
stability  zone  acts  as  a  cap,  preventing  up¬ 
ward  migration  of  methane  and  restricting 
hydrate  formation  to  the  base  of  the  hydrate 
stability  zone,  the  faults  observed  in  these 
data  provide  a  mechanism  in  this  region  for 
migration  of  methane  gas  above  the  BSR. 
Methane  generated  deeper  within  the  thick, 
organic-rich  Blake  Ridge  sediment  may  mi¬ 
grate  upward  through  the  sediment  until  it 
enters  the  hydrate  stability  zone,  where  it 
forms  methane  hydrate,  inhibiting  methane 
from  migrating  farther  upward  and  causing  it 
to  pool  at  the  base  of  the  hydrate-bearing 
layer.  As  sediment  continues  to  be  depos¬ 
ited  on  the  ridge  flanks,  the  added  weight 
causes  slope  failure,  which  is  accommo¬ 
dated  within  the  overlying  rigid,  hydrate¬ 
bearing  sediment  along  the  normal  faults  re¬ 
solved  by  the  seismic  data.  Rupture  along 
these  faults  causes  the  dissociation  of  the 
methane  hydrate  within  the  fault  zones,  en¬ 
abling  significant  amounts  of  methane  to  mi¬ 
grate  out  of  the  gas-rich  layer  and  to  form 
hydrate  within  an  interval  extending  several 
hundred  metres  above  the  BSR. 
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